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We study initial transient stages in directional solidification by means of a non- variational 
phase field model with fluctuations. This model applies for the symmetric solidification 
of dilute binary solutions and does not invoke fluctuation-dissipation theorem to ac- 
count for the fluctuation statistics. We devote our attention to the transient regime 
during which concentration gradients are building up and fluctuations act to destabilize 
the interface. To this end, we calculate both the temporally dependent growth rate of 
each mode and the power spectrum of the interface evolving under the effect of fluctua- 
tions. Quantitative agreement is found when comparing the phase-field simulations with 
theoretical predictions. 

Keywords: Phase-field models; fluctuations; solidification; transient stages. 

1. Introduction 

Solidification processes constitute a primary subject of research in materials science. 
One main objective is to achieve the prediction and control of the final macroscopic 
properties (mechanical, electrical, etc) of solids grown from their melts. These 
properties largely depend on the microstructure, which in turn is result of the con- 
ditions at which the solid was grown [1]. For example, relevant features at the 
microscopic scale such as compositional inhomogeneities or the presence of grain 
boundaries are direct consequences of instabilities and subsequent dynamics of the 
solidification front. Of particular interest for the metallurgical community has long 
been the directional solidification of mixtures and alloys, of direct relevance in pro- 
cesses of zone melting and for the Bridgman method [2] . This configuration is also 
an archetypical model system in nonlinear physics to study pattern selection and 
complex dynamics [3]. In this context it is known that the final wavelength of the 
dendritic array depends on the history [4] , and specifically on the initial destabiliza- 
tion of the solidification front [5,6]. This is a non-steady noise amplification process 
occurring during the early solute redistribution transients, for which considering 
both fluctuations and transient effects are of key importance. Experimental work 
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in this regime can be found in [7, 8] . We are interested here in dealing with these 
aspects in numerical simulations of a directional solidification experiment, by using 
a phase field model supplemented with fluctuations. 

During the last decade, the phase-field model has emerged as a quantitative 
simulation technique to study complex interfacial morphologies [9]. Recently, the 
phase field approach has deserved much attention because it can easily incorporate 
effects like system anisotropies, kinetic attachment or equilibrium fluctuations which 
might be difficult to include with other simulation techniques. This method avoids 
the tracking of the moving front by introducing a continuous order parameter </>(r , t) 
(the phase field), which takes different constant values at the bulk phases separated 
by a smooth interface of width W. The model then consists in a set of coupled 
equations for <fi(r, t) and for the diffusion field which drives the interfacial dynamics, 
and that are constructed so that it reproduces the physical dynamics in the limit 
W — > 0. Early formulations of phase field models were variational, i.e. the model 
equations were derived from functional derivatives of a single free energy functional, 
but the convergence of such models was rather poor. Recent advances in phase field 
formulations, either based in the so called thin interface limit [fO, f f] or in higher 
order expansions [12], permitted to improve dramatically the convergence of the 
model, while dealing with some specific realistic situations. 

However the standard procedure to introduce fluctuations in phase field models 
is restricted to variational formulations [13-15], analogous to the model C of crit- 
ical dynamics [16-18], where the intensity of the fluctuations can be determined 
by using a fluctuation-dissipation relation. As the more recent phase field formula- 
tions do not maintain this variational structure, the fluctuation-dissipation relation 
cannot be used to infer the statistics of the noise appearing in the equations. In 
the model and simulations presented here, we rely on a recent calculus [19] that 
projects the dynamics of a generic stochastic phase field model to the motion of the 
fluctuating interface. This procedure provides a prescription for the intensity of the 
noise terms in the model, accounting for fluctuations of both internal and external 
origin. This approach has previously been used in equilibrium situations only [19], 
and our aim here is to employ it in a transient, out of equilibrium, situation. We 
will restrict ourselves to the study of the initial destabilization in the linear regime, 
where comparison with theory is possible, as a first step to the complete problem 
of considering the whole competition process in the presence of fluctuations. The 
validation of this approach should be also useful in other situations of interfacial dy- 
namics with fluctuations, both in non variational formulations of phase field models 
and in considering more general sources of fluctuations, of external origin. 

2. Stochastic Phase-Field Model for Symmetric Directional Solidifica- 
tion 

In a directional solidification experiment, the mixture sample is pulled at a velocity 
v p in an externally imposed temperature gradient given by T{z) = Tm + Gz, where 
Tm is the melting temperature of the mixture and the tilde refers to physical units. 
For pulling velocities higher than a certain critical value, the interface becomes 
unstable [20] giving rise to cellular regimes. The diffusive field which drives the 
interface dynamics is the reduced solute concentration field u = c ^°° , where Coo 
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is the solute concentration of the sample in the liquid bulk far from the solid-liquid 
interface and ACo = C\ — C% is the miscibility gap which is approximated to be 
constant in this model. We are interested in the symmetric case for which the 
diffusion coefficients in both phases are equal, situation which is typical in dilute 
mixtures of liquid crystals. We consider a modified version of a non-variational 
phase-field model presented by Karma [11,21] for the directional solidification of a 
symmetric dilute alloy. Scaling space and time with I = D/v pi 7 = l 2 /D, being D 
the solutal diffusivity, we obtain the phase-field equations 

ae 2 d t cj ) = £ ^ 2 cj>- f'{cj>)-s\g'mu + ^— -) + eh 1 (v,t) (1) 

l T 

d t u = V 2 u + ^d t h(cj>) - V ■ q(r, t) , (2) 

with 

( V (r,t) V (r' 7 t')) = 2alS(r-r')S(t-t') (3) 

Mr, t')) = 2a 2 u M(r " " *') ■ ( 4 ) 

In the last equations, e = W/l stands for the dimensionless thickness of the interface, 
and It = mACo/G is the thermal length imposed by the external gradient, being 
to the liquidus slope of the alloy. Fluctuations have been included following the 
procedure introduced in [19] to account for the equilibrium statistics. It can be 
shown that the model defined by Eqs. ijTJl. |2"|l. iffil and J3J recovers the classical 
moving boundary description of the problem in the limit of e — > (small interface 
width) and er0,er n <C 1 (small noise intensities). For this particular model, and 
performing a matching asymptotic procedure up to first order in e, we find the 

relations A = a = ^f, a% = ^c^-^n^ > and a l = n°ac$i* ' where J i» 

I2 and I4 are integral constants given by I\ = 2^2/3, I2 — 16/15 and I4 — 0.55, 
and we have chosen the particular set of functions /(</>) = — \(j) 2 + j<p 4 , g{4>) = 
<j) — |0 3 + and h((p) — 4>. In the last relations, (3 is the kinetic attachment, k 
the segregation coefficient, Na the Avogadro's number, v m the molar volume and 
da = gTm I LmACo the chemical capillary length, where a is the surface energy and 
L the latent heat per unit volume. 

3. Simulation Results 

As a quantitative benchmark of the results presented in [19], we perform numerical 
simulations of the phase-field model and compare them with theoretical predictions 
obtained from the moving boundary description of the front. We will present results 
for both the stationary power spectrum of the interfacial fluctuations and for the 
wavelength selection during the initial recoil transient. 

The substance parameters correspond to the liquid crystal 4-n-octylcyanobiphenil 
(8CB), which has been used in several quantitative experimental works [7,22,23]. 
The substance parameters have been obtained from [7,22,23], and are given by 
k = 0.9, to = 88.46 K/mol, T M = 313.5 K, a = 2.2 x 10~ 4 J/m 2 , (3 = 113.04 
s mol/m, D = 4 x 10~ 10 m 2 /s, L = 2.2 x 10 6 J/m 3 , Z = 291.44 g (molecular 
weight) and p — 10 6 g/m 3 (density). The experimental parameters chosen in the 
simulations are = 0.012 mol, v p — 6 x 10~ 5 to/s and G = 2.3 x 10 3 ^C/to. 
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Fig 1. Mean-square interfacial fluctuation spectrum obtained from phase-field simulations and 
compared with the theoretical prediction. 



3.1. Stationary power spectrum 

We have simulated a planar, stationary interface in order to evaluate its spec- 
tral statistical properties. In particular, we consider the mean-square fluctuation 
spectrum, which can be predicted analytically by S(k) = f ^-(zf(k)zf(k')) = 
NA{i-k)(G+md k 2 ) ' being Zf(k, t) = J dk Zf(r, t) e _lkr the Fourier transform of the 
interface position. Figure ^ compares the theoretical power spectrum S(k) with 
the power spectrum obtained by analyzing the spectral properties of the phase-field 
interface. This simulation has been performed by using an explicit finite differ- 
ences scheme in a N x x N z = 256 x 100 rectangular grid with a scaled time step of 
At = 6.75 x 10~ 5 and a scaled grid spacing Az = Ax = 0.03. The interface width 
is e = 0.0375, and the phase-field simulations have been thermalized during 10 4 
time steps before calculating any of the magnitudes presented in this work. The 
interfacial fluctuation spectrum has been averaged among 10 samples taken every 
10 simulation steps in order to avoid statistical time correlations. The quantitative 
agreement indicates the correctness of the determination of the model parameters. 
The scale associated to the interfacial width e, which constitutes the limit of the 
descriptions of any phase-field model, is indicated by a vertical dashed line and 
determines the maximum wavelength for which the simulations results are reliable. 

3.2. Transient power spectrum and wavelength selection 

We next investigate the amplification of the fluctuations during the initial redis- 
tribution transient. The simulations are compared with a theoretical prediction 
based on the Warren and Langer's work on noise amplification [5]. In this ref- 
erence, the time evolution of the interface correlations is described by S(k, t) = 

So(k) J^oq cte'exp |2 J* t * Lu(k,t)dt'^ where Lu(k,t) is the transient growth rate of the 

symmetric model [24], and So(fc) is given by S'o(fc) = lgIi^^Jt fc 2 ) ' being uj e — 
to(k,t = oo) the amplification rate at t — > oo. Figure [21 shows the amplification of 
the interfacial fluctuation spectrum at three different physical times during the ini- 
tial transient (f= 0, 1.8 s and 2.7 s). Dashed lines are obtained by integrating S(k, t) 
in the adiabatic regime u> <C fc 2 . In this case the interface thickness is e — 0.15, and 
the simulations have been performed in a 500 x 256 grid with Az — Ax = 0.09, 
At = 0.00135. The resulting power spectrum is noisy because the results are di- 
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Fig 2. Noise amplification and wavelength selection during the transient. Comparison between 
the theoretical prediction and phase-field simulation results. 



rect spectral properties of the fluctuating interface at each time and no average or 
filtering has been used. As it can be observed, quantitative agreement is achieved 
between theory and simulations in the early time amplification of fluctuations as 
well as in the wavelength selection of the cellular pattern. After the early stages 
where the growth of fluctuations is linear, the system enters in a nonlinear regime 
and a cellular pattern with a certain wavelength sets in the system. The end of the 
linear regime can be defined as the time at which the amplitude of the most unstable 
Fourier mode is comparable with its wavelength. This condition defines a crossover 
time to which can be theoretically determined by the condition < Sz 2 >~ X ma x, 
where < Sz 2 >= f S(k, t) is the mean-square fluctuation amplitude and \ max 
is the wavelength of the largest Fourier mode. Figure shows the evolution in 
time of these two magnitudes in our system, and indicates the determination of the 
crossover time to i n physical units. 




Fig 3. Determination of the crossover time to. 



4. Conclusions 

We have studied the initial transients of a directional solidification experiment of 
a mixture of liquid crystals in the context of the symmetric model with constant 
concentration gap. Contrary to previous work on transients in fluctuating phase- 
field simulations [24], results presented here have been obtained by using a non- 
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variational formulation of the stochastic phase-field model, recently proposed [19]. 
Fluctuations have been implemented giving the correct equilibrium statistics but 
without employing any fluctuation-dissipation relation to determine the noise in- 
tensities in the model. The agreement between predictions and simulations has 
been excellent, thus confirming the applicability of this approach for the study of 
fluctuations during transient stages. As a final conclusion, the employ of this class 
of efficient stochastic phase-field models opens new possibilities for the quantitative 
study by simulation of microstructure formation in solidification processes. 
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